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ABSTRACT 

The galaxy cluster RX J1347.5— 1145 is one of the most X-ray luminous and most massive clusters 
known. Its extreme mass makes it a prime target for studying issues addressing cluster formation and 
cosmology. Despite the naive expectation that mass estimation for this cluster should be straightfor- 
ward (high mass and favorable redshift make it an efficient lens, and in addition it is bright in X-rays 
and appears to be in a fairly relaxed state), some studies have reported very discrepant mass esti- 
mates from X-ray, dynamical and gravitational lensing. In this paper we present new high-resolution 
HST/ACS and Chandra X-ray data. The high resolution and sensitivity of ACS enabled us to detect 
and quantify several new multiply imaged sources, we now use a total of eight for the strong lensing 
analysis. Combining this information with shape measurements of weak lensing sources in the central 
regions of the cluster, we derive a high-resolution, absolutely-calibrated mass map. This map pro- 
vides the best available quantification of the total mass of the central part of the cluster to date. We 
compare the reconstructed mass with that inferred from the new Chandra X-ray data, and conclude 
that both mass estimates agree extremely well in the observed region, namely within 400/ifo^kpc of 
the cluster center. In addition we study the major baryonic components (gas and stars) and hence 
derive the dark matter distribution in the center of the cluster. We find that the dark matter and 
baryons are both centered on the BCG within the uncertainties (alignment is better than < 10 kpc). 
We measure the corresponding 1-D profiles and find that dark matter distribution is consistent with 
both NEW and cored profiles, indicating that a more extended radial analysis is needed to pinpoint 
the concentration parameter, and hence the inner slope of the dark matter profile. 
Subject headings: cosmology: dark matter - gravitational lensing - galax- 

ies:clusters:individual:RX J1347.5-1145 



1. INTRODUCTION 

Galaxy clusters have been the focus of very intense re- 
search over the past decade. They are extremely valuable 
for studying the empirical properties of dark matter, as 
well as for exploring the growth of structure at the high 
mass tail of the mass function. The mass distribution of 
galaxy clusters is particularly important for cosmologi- 
cal studies because it provides a critical test of the Cold 
Dark Matter (CDM) paradigm and constrain dark en- 
ergy models, if their mass and its distribution is reliably 
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determined. However, in some cases the cluster mass es- 
timates from the weak lensing, strong lensing and X-ray 
measurements disagree: the origin of these discrepancies 
is likely to be a combination of erroneous assumption of 
the cluster potential, projection effects and complicated 
gas physics not taken into account. In addition, sys- 
tematic effects like cluster member contamination and 
unknown redshift distribution of sources, strong lensing 
image identification, temperature calibration, etc. enter 
the error budget for X-ra y and lensing analyse s. 

It was first proposed bv I Navarro et al.l (|l997l ) that the 
dark matter halos on a variety of scales should follow 
a universal profile within the currently accepted ACDM 
paradigm. The 3-D density distribution of dark matter 
should follow pdm oc within a scale radius and falls 
of steeper at radii beyond that (pdm cx r~^). However 
recent simulations have shown that the profile might be 
modified and th e structural paranieters vary with red - 
shift and mass fNeto et al.l l2QQ7l : iNavar ro et all l2QQ4f ). 
In addition, likely interaction between da rk matter and 
baryons further complicate t he picture (|Gnedin et al.l 
I2QQ4I : iNagai fc Kravtsovl[2QQ5h . It is therefore crucial to 
compare simulations and observations keeping the bary- 
onic component in mind. This will allow us to mea- 
sure the slopes of dark matter profiles, which is a crit- 
ical test for our currently accepting cosmology as well 
as understanding the complicated baryonic physics in 
galaxy clusters. Several works have previously stud- 
ied mass distribution in number of clusters using com- 
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Nataraian & Kneiblll996l:lKneib et alJl2003l: ISmith et al. 
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Dieffo et all l2005l: ICacciato et alJ 120051: iJee et al. 


20071; 


Limousin et al.l |2007f) and combined strong lens- 




ing and stellar kinematics data of the dominating cen- 
tral galaxy ( S and et al. 2007). These approaches offer 
valuable extra constraints for determining the mass dis- 
tributions. 

In this paper we study the most luminous X-ray clus- 
ters known to date. The cluster RX J1347.5-1145 
(Schindler et al. 1995) at a redshift z = 0.451 has been 
a subject of intense past research. Rich datasets includ- 
ing X-ray (iSc hindler et al."1995', ^997"; "Allen et al."2002'; 
Ettori et alT T2004: Gitti & Schindler 2004; Gitti et al. 
and optical (iF ischer & Tvson 1997; Sahu et al. 
iCohen fc Kneib 2002; Ravindranath fc Ho 2002) 
as we ll as observations of Suny a ev-Zerdovich (SZ) ef- 
fect (iPointecouteau et aH l2001l : iKomatsu et"al] l2001l : 
iKitavama et al.l 120041 ) have been obtained, yet the mass 
determinations based on X-ray properties, SZ effect 
(SZE), velocity dispersion measurement, strong and 
wea k lensing have many tim es yielded discrepant results 
(see I Cohen fc KneibI l2002l for a summary). In partic- 
ular th e discrepancy between the dynamical mass es- 
timate (Cohen fc Kneib" "2002), early X -ray mass mea- 
surements (Schindler et al. 1997) and lensing results 
([Fischer fc Tvsonlll997l ) yielded a factor of ~ 3 discrep- 
a ncy in their mass estimates. 

ICohen fc KneibI (|2002l ) suggested the cluster is likely 
undergoing a major merger which would reconcile the low 
velocity dispersion (and hence virial mass) measurement 
and its high mass as predicted by gravitational lensing 
and X-ray data. If the velocity dispersion measures pre- 
dominately one component, and the X-ray gas has viri- 
alized to the post-merger cluster mass, this can possibly 
explain the data. Further evidence for a merger scenario 
comes also from the optical morphology: the cluster con- 
tains two cD galaxies (see e.g. Fig. [T] - throughout this 
paper we will refer to the brighter, western one as the 
BCG, and the other as the second cD). Finally, there is 
a region of shocked gas in the so uth-ea st quadrant dis- 
covered and d iscussed by iKomatsu et al.l ()2001l ) (with 
SZ effect) and lAllen etlo] (|2002h based on X-ray obser- 
vations with Chandra. We confirm the latter using the 
newest Chandra X-ray data. 

In order to r esolve the puzzle of mass estimates 
in IBradac et al.l (|2005al ) we studied this cluster using 
ground-based optical data, and for the first time de- 
termined its mass distribution using a combined strong 
(multiply-imaged systems) and weak (statistical mea- 
surement using shape information of an ensemble of 
background galaxies) gravitational lensing analysis. The 
main drawback of that reconstruction was the lack of 
clearly identified multiply imaged systems and their red- 
shifts. In the current paper we obtained new multi-color 
HST/ACS data. The high resolution and sensitivity of 
ACS allowed us to unambiguously identify many new 
multiply imaged systems. In this analysis we are able to 
use 9 multiply imaged syst ems, in contrast to the single 
one of lBradac et al.l (|2005af ). Furthermore the far greater 
density of sources that can be used for weak lensing (a 
factor of 5 improvement) has allowed us to obtain an 
absolutely-calibrated mass map at high resolution. We 
also present new X-ray data and compare estimates from 



both analyses. In addition we describe a new automated 
search routine for multiply imaged sources we are cur- 
rently developing; its main advantage is that it combines 
colors and information on lensing geometry of the whole 
system to obtain matching pairs in a semi-automated 
way. 

Our final mass reconstruction is again performed using 
a pixelated model which combines the information from 
strong and weak gravitational lensing. If spectroscopic 
redshift s of the multiply imaged systems were available, 
this would fully eliminate the need of using models that 
assume a particular shape of the gravitational potential. 
However, since only one system has a spectroscopically 
confirmed redshift, we use the predictive power of gravi- 
tational lensing to predict the redshifts for the remaining 
systems. As explained further in ^ we therefore model 
strong lensing data using parametrised models at the 
start to obtain the redshifts. For the final reconstruction 
we however do not assume a special shape of the gravi- 
tational potential, instead exploring a much larger range 
of cluster mass distributions and profiles (provided the 
estimated redshifts are reasonably reconstructed). Fur- 
thermore, we use cluster to study the interplay between 
baryons and dark matter, therefore we study separately 
its stellar, gas and dark matter component. This allows 
us to measure the cluster's dark matter profile, one of 
the strong predictions of ACDM cosmology. 

This paper is structured as follows. In section [2] we 
describe the optical data used in this analysis. The basic 
image processing and the extraction of the strong gravi- 
tational lensing data is described in section [3] and in sec- 
tion m we describe the weak gravitational lensing data. 
We present the X-ray data in section[5l We infer the mass 
distribution of the cluster RX J1347.5— 1145 from lensing 
data, and compare it with X-ray data, in section [6l The 
individual contributions of stars, gas, and dark matter 
to the total mass of the cluster are studied in section [71 
The conclusions are presented in section [51 To evaluate 
the angular diameter distances throughout the paper we 
assume the ACDM cosmology with = 0.3, I^a = 0.7, 
and Hubble constant Hq = 70 kms~^ Mpc~^. 

2. OPTICAL OBSERVATIONS AND INITIAL DATA 
REDUCTION PROCESS 

ACS/WFC imaging of the cluster RX J1347.5-1145 
was carried out in Cycle 14 (proposal 10492, PI Erben) 
on 2006 March 9-11. The cluster was observed in a sin- 
gle (dithered) pointing and three different filters F475W, 
F814W, and F850LP for 5280s (two orbits) each. The de- 
mands placed by the lensing analysis require special care 
when reducing the images. We base the reduction on the 
bias and flat-field corrected fit-images (provided by the 
standard ACS pipeline). For each exposure we compute 
a noise model including all noise sources except object 
photon noise, and update the bad pixel mask as detailed 
in Marshall et al. (in preparation). We subtract the 
sky background seperately in the four image quadrants 
to correct for the sometimes present residual bias level. 
Satellite trails are masked manually, in order to exclude 
them from the co addition. We use the Multidrizzle 
([Koekemoer et al.| [2002) routine to align the images, cor- 
rect for geometric camera distortion, mask cosmic rays, 
and for coaddition. To register the images with the as- 
trometric accuracy needed for the lensing analysis, we 
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determine the offsets among the images by m atching win- 
dowed SExtractor (|Bertin fc Arnoutsll 19961 ) positions of 
high S/N objects in the individual, distortion-corrected 
exposures. We calculate residual shifts and rotations us- 
ing the IRAF routine geomap, which are then fed back 
into Multidrizzle. We use "square" as the final driz- 
zling kernel, where we slightly shrink the input pixels 
(pixfrac = 0.9) and set the output pixel scale to 0.03 
arcsec. This is smaller than the original pixel scale of 
the ACS/WFC detector in order to reduce the impact 
of resampling on the shape measurements. In the final 
coaddition, we use our noise model for inverse variance 
weighting. For further details on the data reduction see 
Marshall et al. (in preparation). 

3. STRONG LENSING DATA ANALYSIS AND FITTING 
PARAMETRISED MODELS 

With the advent of the high resolution HST optical 
images the quality and complexity of the strong lens- 
ing data increased dramatically. Many clusters have 
been observed with HST showing numerous multiply im- 
aged systems (the most notable exam ple b eing A1689 - 
iBroadhurst e t alll2005l : UTalkola et"an(200^ . By far the 
most time-consuming part of the strong analysis is, how- 
ever, to match the multiply imaged system - i.e. identify 
all the multiple images that were produced for a single 
source. 

First steps in trying to op timise this process w ere done 
bv iKneib etHI (p.993) and iSharon et all (|2005f ). Their 
methods still involve significant human involvement; we 
are aiding the process of searching for multiple image sys- 
tems with automated color and geometry matching. Our 
method is currently under development and once finished 
will allow for simultaneous search of the best matching 
objects in the color space as well as including the ge- 
ometry of the system (i.e. positions and orientations of 
the lensed images). This is a very important step in the 
view of large samples of galaxy clusters that will be dis- 
covered and targeted in the future optical surveys (such 
as SNAP, LSST, DUNE and PanSTARRS). 

We start by considering a simple initial guess for a 
cluster potential together with a hypothesis of at least 
one multiple image system. We use the image positions 
to determine the mass distribution of the cluster using 
a simple parametric model. The image positions (at 
the end of this process we use all images listed in Ta- 
ble □ and shown in Fig. [1]) are used in LENSTOOL, 
(http : // www . oamp . f r/ cosmology/lenstool/, 
iKneib et al. 1993). The code determines the best 
fit parameters of a parametrised model. In addition, 
its Bayes optimization and Monte Carlo Markov Chain 
(MCMC) sampling routines (Julio et al. 2007) allow us 
to obtain a set of all models that satisfactory fit the data 
(for few model parameters, such as the mass and ellip- 
ticity of the cluster). We characte rize RX J1347.5— 1145 
by two PIEMD profil es (see e.g. iLimousin et al.l [20051 : 
lElfasdottir et al.l 12007 for the basic properties of the 
model), centered on the two brightest cluster members. 
In addition we also include the 20 brightest cluster 
members in the I-band to the mass model. They are 
each modelled as PIEMD spheres with a line-of-sight 
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scaling was found to be ctpiemd ~ 

(260 ±50) kms"^ and 
^cut PIEMD = (5 ± l)kpc for a galaxy with mj = 20.5. In 
total this adds to 50 constraints and 34 free parameters 
(42 when including unknown redshifts). 

The parameters of the best-fit model (using the final 
catalogue of all multiply imaged systems) are listed in 
Table [2l The ellipticities of the two main components 
are given as the ratio of minor/major axis {b/a) and 
the position angle (j) measured west-through-north. The 
parametrised model predicts the images with an average 
accuracy of 3''. The least well predicted image system is 
the image system F, which is however perturbed by one 
of the cluster members (included in the model only as 
one of the 20 cluster members and not minimized for in- 
dividually) . The two main components have very similar 
velocity dispe r sions, however as discussed extensively in 
IHalkola et al.l (|2008l ) given the degeneracies in the mod- 
elling, this is not the only solution for this cluster. 

We justify the choice of a two component model 
by comparing the Bayesian evide nce (for details see 
IKneib et al.ll2003l : lJullo et al.ll2007[ ) when fitting one or 
two PIEMD components (both times including the clus- 
ter members to the fit, and letting all the parameters 
of one or both PIEMD components to vary). The log- 
arithmic evidence is a factor of 3 larger in the case of 
a two component model, thereby clearly justifying the 
two-clump model used for the fit. 

Using the samples (i.e. a set of models that satisfactory 
fits the data) we then predict possible additional images 
for the initial guess system and matches to further (not 
necessarily yet identified as multiple) images. This allows 
us to find a region in space, where matches to known 
systems and objects should be searched for. Further, if 
a match in color space is found, the hypothesis that it is 
a multiply imaged system can be tested and the redshift 
estimate is obtained using the improved lens model. 

All the identified images were visually inspected (we 
are hoping to be able to omit this step altogether in the 
future). 

3.1. Redshifts for Strongly Lensed Systems 

For RX J1347.5— 1145 the only multiply imaged system 
with confirmed spectroscopic redshift is system A (Ta- 
ble [1] and Fig.[T]). The redshift of the brighter of the two 
images was independently measured to be 1.75 for both 
images with FORS spectroscopic data (|Lombardi et al.l 
I2OO7I ) as well as Keck data. This is in agreement with 
the photometric red shift obtained from UBVRIJHKs 
ground-based data in lBradac et al.l (|2005al ). where a red- 
shift 1.76 ± 0.1 was used for this system. We have pre- 
dicted the redshifts of the other identified systems by set- 
ting them as free parameters when the system was added 
in the image identification procedure outlined above. 
The resulting best-fit redshifts are given in Table [TJ The 
final error bars were obtained by source-plane modelling 
(modelling in the image-plane with all the parameters set 
free would be too time consuming) and should be treated 
as approximate. 

Ideally, the unknown redshift should be fitted simul- 
taneously with weak and strong lensing data using the 
method described in ^ (which does not assume a spe- 
cific model family for the cluster potential). This is un- 
fortunately not possible due to the so-called mass-shee t 
degeneracy (jFalco et al.l 119851 : ISchneider Seitd Il995f ). 
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TABLE 1 
The properties of the 
multiply-imaged systems used in this 
WORK. The redshift of the system A 

WAS OBTAINED USING FORS SPECTRA, 
WHEREAS THE OTHERS HAVE BEEN 
DETERMINED BY LEAVING THEM AS FREE 
PARAMETERS WHEN FITTING STRONG 
LENSED DATA WITH PARAMETRISED 
MODEL. 





Ra 


Dec 


^pred 




zUo.oYzUY 


— ll.YblOTz 




A 


zUo.oozoo 


— ii. m44U 1 


1.75 




206.87250 


-11.746470 






206.87117 


-11.748356 




B 


206.87316 


-11.745810 


1.2 ± 0.1 




206.87423 


-11.745110 






206.88293 


-11.741210 




C 


206.88415 


-11.741798 


2.0 ±1.0 




206.87845 


-11.749200 




D 


206.87833 


-11.749682 


2.2 ±0.1 




206.87219 


-11.765260 




E 


206.87406 


-11.766258 


2.5 ±0.6 




206.86549 


-11.764203 




F 


206.86661 


-11.765704 


4.0 ±2.0 




206.86395 


-11.755330 






206.87357 


-11.768255 




G 


206.87130 


-11.767085 


3.0 ±0.5 




206.88088 


-11.770223 




H 


206.86634 


-11.766999 


4.2 ±1.0 




206.86714 


-11.742389 






206.88754 


-11.757578 






206.88513 


-11.748397 




I 


206.87878 


-11.753372 


1.7 ±0.2 




206.87765 


-11.759400 






206.86960 


-11.747372 







This transformation essentially allows us to determine 
the surface mass density only up to a constant (that en- 
ters in a re-scaling and offset of the surface mass density). 
Since the redshift information enters the len sings; equation 
in similar fashion (see Br adac et alJ |2QQ3 for more de- 
tails), if only a single redshift of multiply imaged system 
is known, without assuming a parametrized model, the 
rest of the redshifts can not be determined. Weak lens- 
ing in principle should provide a second redshift plane 
and hence constrain the profile and its normalisation, al- 
lowing us to reconstruct the data without the need of 
any simple parametrised modelling. In practice, how- 
ever, the weak lensing data is noisy and therefore not 
sufficient to break the mass sheet deg eneracy (see e.g. 
iBradac et aIll2QQ4l : iLimousin eraII[2QQ7 ). If one tries to 
constrain the unknown redshift of a system with images 
far away from the images with known redshift, one can 
obtain equally good fits by either slightly perturbing the 
potential or by changing the redshift. Hence we need 
to use parametrized models for this part of the analy- 
sis. However, when other multiply imaged systems are 
present in the vicinity, the redshift can still be deter- 
mined using the pixelated model; when we added image 
system I we could only get a good fit with z = 1.7 and 
not the photometric redshift estimate of z = 2.2; the 
same conclusion was obtained already by iHalkola et al.l 
[20M 

4. WEAK LENSING DATA ANALYSIS 

We perform the weak lensing shape measurements 
on the F814W data due to their higher galaxy num- 



TABLE 2 

Results of best-fit parametric modelling of strong lensing 

DATA (parameters REFER TO THE VALUES USED BY LENSTOOL, 
AXIS-RATIO b/a AS THAT OF THE PROJECTED MASS DISTRIBUTION). 



(^PIEMD Tc Vcut b/a (p 

( kms~^) (kpc) (kpc) 



BCG 1640 ± 40 340 ib 30 1000 ib 30 0.89 ± 0.01 16° ± 1° 
2nd cD 1640 ± 10 400 ib 5 2000 ib 60 0.87 ± 0.01 70° ± 1° 



ber density. For the PSF anisotropy and smearing 
correction we clo s ely fo llow the technique described in 
ISchrabback et al.l (l2QQ7f) The p rocedur e is based on the 
KSB + algorithm (Kaiser et al. 1995l; iLuppino fc Kaiserl 
Il997l : fHoekstra et al . 1998 ), which has extensively been 
tested on the simulations of the STEP collaboration 
(|Hevmans et al.ll2QQ5l : iMassev et al.ll2QQ7[ ). The KSB+ 
method is formally valid only in the weak lens ing regime, 
however as discussed i njBradac et al.l (|2QQ 6) (see also 
iMa ssev & Goldberg 200^ the approach is also valid in 
the non-weak lensing regime to the accuracy needed here. 
We select stars from the half-light-radius vs. magnitude 
diagram. 

In contrast to Bradac et al. (2006) we use the new PSF 
correction scheme developed by Schrabba ck et al. ] ( |2QQ7| ). 
The method is superior, since it takes the temporal (as 
well as the spatial) variability of the ACS PSF into ac- 
count. The PSF variation is measured from stars in indi- 
vidual exposures (rather than in the co-added image as 
was traditionally done). It is then compared with star- 
rich fields from archival observations (stellar templates). 
The shapes of the stars in individual exposures of our 
data set are first compared to the shapes of stars from 
stellar templates and the best template is chosen. The 
full PSF model is then composed of several templates; 
such a model therefore captures the full spatial as well as 
temporal variability of ACS. We detect short-term vari- 
ations of the ACS PSF, which are interpreted as focus 
changes due to thermal breathing of the telescope (see 
also Schrabba ck et al. 2007; Krist 2003; Rhodes et^ 
l2007l : lAnderson fc Kindl2006f ). The detailed PSF model 
obtained from this procedure is then used to correct the 
shapes of all the objects in the RX J1347.5— 1145 field. 

Finally, since the PSF shape of stars changes as a func- 
tion of the scale at which we use for the sha pe measure- 
ment (see e.g. iHevmans et al.l 120051 : iJee et al.l 120061) . we 
match the scale used for PSF and galaxy shape mea- 
surement. This method reduces the systematic contri- 
bution to the shear correlation functions due to PSF 
distortions to < 2 x 10~^ for galaxy fields containing 
at least 10 stars, and was demonstrate d to be sufficient 
even for cosmic shear measurements (Schrabback et aP 
l2007l ). where the signal is more than an order of magni- 
tude smaller than measured in cluster fields. In addition, 
we apply a parametric model, similarly to the one pro- 
posed by Rhodes et al. (2007), to correct for the infiu- 
ence of the degraded charge-transfer-efficiency on galaxy 
shapes, see Schrabback et al. (in prep.) for further de- 
tails. Us in^ simulated weak len sing data from the STEPl 
project (|Hevmans et al.l [20061 ). we identified a constant 
bias between the input and measured shear values. In 
the analysis of the STEP2 data we correct for this bias 
by the introduction of a multiplicative shear calibration 
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Fig. 1. — The F450W-F814W-F850LP color composite of the cluster RX J1347.5-1145. Multiply imaged systems are marked and labeled 
(see also Table [T}. North is up and East is left, the field is 2.3^ x 2.3^, which corresponds to 800 x 800 kpc at the redshift of the cluster. 
The color composite was created following the algorithm from iLupton et al.i (^2004i) . 



factor of 1.10, which prove to be on average accurate 
at the ^ 2 % level also for thi s different set of image 
simulations (Massev et al J [20071 ). In the analysis of the 
RX J 1347. 5— 1145 data we apply slightly different galaxy 
selection criteria, leading to a marginally changed cali- 
bration factor of 1.08, which we apply to the data. 

To exclude cluster members from the final catalog 
we match it to the photometric redshift catalog from 
the ground-b a sed (V LT/FORS and ISAAC) data from 
iBradac et al.l (|2005a| ). We exclude all the galaxies hav- 
ing measured photometric redshifts ^phot < 0-5. Un- 
fortunately, however, we can not use the photometric 
redshift estimates for all the object in the weak lensing 
catalog, since many of the faint objects were undetected 
or unresolved in the groun d-based data. Th e refore fol- 
lowing the prescription of ISchrabback et all (|2007f ) we 
use the GOODS-MUSIC sample and apply the corre- 
sponding magnitude cuts to the data and obtain the aver- 
age redshift for galaxies in the weak lensing catalogue of 
zwL = 1-4. The exact value we use has little importance 
(as discussed in §6.3p on the final mass reconstruction, 
since the weak lensing signal/noise is much lower than 
the strong lensing one, and therefore the mass scaling is 
determined predominately by the strong lensing systems. 
This is however only true for the regions where we recon- 
struct the mass, once we extend the measurements far 
outside the strong lensing regime with the ground based 
data (beyond ~ 500 kpc) the determination of redshift 
distribution will become much more important. 

5. X-RAY DATA REDUCTION 

Under the assumptions of hydrostatic equilibrium and 
spherical symmetry, one can use the observed X-ray sur- 
face brightness profile and the deprojected X-ray gas 
temperature profile to determine the total mass and gas 



mass profile. The reduction and anal ysis of the Chandr a 
X-ray data is described in detail bv I Allen" et all (|2007r ). 
ISchmidt fc AllenI (|2007f ). and Million et al. (in prepa- 
ration). The thermodynamic X-ray pressure map de- 
termined from the Chandra data (total clean exposure 
time 67.9ks) is shown in Fig. [3l Overall, the pressure 
map shows that the cluster is remarkably relaxed, al- 
though the localized region of high-pressure, shocked 
gs:as in the southea s t qua drant d iscovere d and discussed 
by iKomatsu et all (|200lh and lAllen et all ((2002) is also 
clearly visible (indicating recent merger activity). 

The cluster mass profile was determined from the ob- 
served deprojected X-ray temperature (Fig. [S]) and sur- 
face brightness profiles (see Schmidt & Allen 2007 for de- 
tails). The data from the southeast quadrant (position 
angles of 180-280 degrees measured west-through- north) 
were excluded from this analysis. In order to compare 
with the lensing mass measurements presented in this pa- 
per, we have projected the three-dimensional mass (and 
later also gas) profile determined from the Chandra data, 
assuming that the cluster extends to the virial radius (as 
determined from the X-ray data), i^vir = 3.03 Mpc. To 
estimate the systematic errors arising from this choice 
of truncation radius, we adopt conservative limits and 
change the radius by 30%. This resulted in changes to 
the projected masses measured within 600kpc of the clus- 
ter center that are within the quoted uncertainties; the 
precise value for the truncation radius is therefore of lim- 
ited importance for this study. The resulting projected 
integrated mass profile is plotted in Fig. [5] together with 
mass estimates obtained from lensing analysis. 



6. CLUSTER MASS RECONSTRUCTION FROM STRONG 
AND WEAK LENSING DATA 
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Fig. 2. — 10^^ X 10^' cutouts of the multiply imaged systems from the 
BRK FORS data was used due to its extreme red color). 

In this section we present the combined strong and 
weak lensing mass reconstruction of RX J1347.5— 1145. 
Here we ess entiahy fohow the method described in 
iBradac et al.l (2QQ5b) and implemented on the ACS data 
in ' Bradac et al.i (^2006|). The basic idea is to describe 
the cluster's projected gravitational potential by a set of 
its values on a regular grid i/jk^ from which we evaluate 
all quantities relevant for gravitational lensing by finite 
differencing. E.g. the scaled surface mass density is 
related to ip via Poisson equation, 2tv = V'^ip (where the 
physical surface mass density is H = Scrit and Scrit is 
a constant that depends upon the angular diameter dis- 
tances between the observer, the lens, and the source). 

The advantage of such an approach is that we avoid 
making any assumption on e.g. shape and/or profile of 
the potential, which is crucial when dealing with merg- 
ing clusters. The strong and weak lensing data are then 
combined in a x^-fashion. We minimise the by search- 
ing for the solution of the equation dx^/dipk = 0- We 
linearize the resulting system (using the initial values for 
the potential - see below - or values from the previous 
iteration) and solve it using sparse matrix techniques. 
Since the weak lensing data is noisy, regularisation (i.e. 
a process which ensures that unphysical pixel-to-pixel 
variations in surface mass density are suppressed) needs 



F450W-F814W-F850LP color composite (except for image F, where 

to be employed. For this purpose we compare the cur- 
rent surface mass density map with one obtained on a 
coarser grid, and pe nalize strong deviations in (see 
iBradac et al .1 l2QQ5bl for details). The regularisation is 
chosen (and tested to be) such as to not bias the result- 
ing mass estimates. 

6.1. Initial conditions of the method 

For the purpose of obtaining the initial values of the po- 
tential needed in the first step of the iteration we assume 
a single singular isothermal sphere centered on the BCG 
with velocity dispersion a = 1600 km s""*^, taken from 
the st rong lens modeling. As noted bv 'Elfasdottir e t al] 
(I2007D there is a non-trivial conversion between a and 
c'"piemd; therefore we simply assume a model with sim- 
ilar enclosed mass within the Einstein radius. We have 
chosen a more simplistic model than in ^ on purpose, 
in order not to bias our results and obtain a reconstruc- 
tion independent of the assumptions about the potential 
(as is the case in assuming parametrised models). These 
particular assumptions, however, do not infiuence the re- 
sults much, when high quality data, such as the ones here, 
are used (see discussion in §6.3p . Therefore any reason- 
able guess for a cluster potential (from X-rays, dynamical 
mass estimate, etc) can be used. 
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Fig. 3. — (top) The deprojected temperature profile of the X- 
ray emitting gas in RXJ1347. 5-1145 determined from the Chandra 
X-ray data. The data from the southeast quadrant were excluded 
from the analysis, as described in the text, (bottom) The pres- 
sure map of RX J1347.5— 1145. The two crosses indicate the two 
brightest cluster members. 

6.2. Combined reconstruction 

As described in section [3] we use 9 different multi- 
ply imaged systems and 700 weakly-lensed galaxies 
(70 arcmin"^). We start with a 30 x 30 pix^ grid for a 
4.2 x4.2 arcmin^ field (8'Ypix; the cluster is not centered 



in the ACS field) oriented with respect to the negative 
RA coordinate (not all of the grid cells in the field con- 
tain data and those that do not are excluded from the 
reconstruction). We gradually increase the number of 
grid points in steps of 1, with the final reconstruction 
performed on a 60 x 60 pix^ (4'Ypix) grid. 

The resulting reconstruction is shown in Fig. [3J to- 
gether with X-ray surface brightness contours from the 
Chandra observation (see Sect.[5|). We clearly detect the 
main cluster component, which is aligned with the BCG. 
The offset from the BCG is {-2" ± 2" , -I" ± ?>") and is 
also consistent with the peak of the X-ray surface bright- 
ness. We also resolve the south west structure, which also 
shows an overdensity of cluster members. The resulting 
model predicts the strongly lensed image positions with 
an average accuracy of less than A" ^ which is also effec- 
tively the final pixel size. In order to reach higher resolu- 
tions, adaptive grid methods needs to be employed; this 
will be a subject of future work. 

The mass profile, calculated by determining the en- 
closed, projected mass within circular apertures is plot- 
ted in Fig. [5l The errors include both systematic and 
statistical contributions, their estimation is described 
in more detail in Sect. 16.31 We fit power-law model 
(M(< R) oc R^) to the profile, and find it to poorly 
represent the data. The logarithmic slope of the best 
fit model is n = 1.6 ± 0.1 (shallower than isothermal). 
Adopting an isothermal profile, we find a line-of-sight 
velocity dispersions for a = (1550 ± 100) km s~^. 

The total mass profile is in excellent agreement with 
the X-ray data . It also agrees with the analysis 
in iBradac et al.l (l2Q05af) . and strong lens modeling in 
iHalkola et al.l ( 20081 ) and the one presented in ^ (the 
data points using only strong lensing information in Fig [5] 
are extrapolated to a larger radius using the adopted pro- 
file for modelling). The agreement between our strong 
and strong and weak lensing modelling is expected, since 
for the analysis i n ^ we use t h e sam e set of multiply 
imaged systems; iHalkola et al.l (|2008f ) does use a sub- 
set of multiple images, and the agreement between these 
semi-independent analyses is very good. Most impor- 
tantly, however, the agreement between X-ray data and 
lensing analyses was not achieved to this accuracy in the 
past and it is very encouraging to note that with ex- 
cellent data we can indeed measure cluster masses re- 
liably. Whereas the results from Bradac et al. (2005al) 
show that even with the ground-based data one can re- 
liably measure enclosed mass, in order to obtain the full 
mass profile space-based data are required. In addition, 
as shown by IHalkola et al.l (|2008l ) , if using a reasonable 
guess for the cluster potential only a few systems are 
needed to reconstruct the mass profile. However, further 
image systems, and the addition of weak lensing data, 
increases the resolution and radial extent of the recon- 
struction, which helps breaking the degeneracies in the 
mass profile as shown in ^ 

We also compare our estimates with the recent X- 
ray an d gravitational lensing comparison of Gitti e t al.l 
(|2007[ ). The paper is based on old lensing data (the 
ACS/HST data is not included) together with XMM 
data. Whereas the lensing results have not changed sig- 
nificantly, there is a disagreement between Chandra and 
XMM X-ray data. 
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The X-ray results presented here are in excellen t agree- 
ment with those reported by lAllen et all (|2QQ2l ) from a 
previous analysis of a Chandra data set with shorter 
exposure. The X-ray derived masses are larger than 
those reported from the analysis of XMM- Newton data 
by Gi tti et al . (2007). The origin of this discrepancy is 
likely to lie at least in part in the complicating effects 
of the XMM-Newton point spread function, which are 
significant in the presence of strong temperature gradi- 
ents (Fig, pi ) and which were not modelled in detail by 
iGitti et al.ll2Q07l 

The largest circular aperture we can place on the BCG 
is 350 kpc. The enclosed, projected mass within this 
aperture is Ml(< 350 kpc) = (5.9 ± 0.5) x IO^^^Mq. The 
uncertainties include statistical as well as systematic er- 
rors, originating from weak lensing measurements, pos- 
sible multiple image and redshift misidentification, and 
the initial model we used. We now describe the error 
budget in more detail. 

6.3. Errors and Possible Systematic Effects 

As in iBradac et al.l (|2006l ) , we also extensively study 
possible errors and systematic effects. We generate 1000 
bootstrap resampled weak lensing catalogs and perform 
reconstructions on each of these. To further test the re- 
liability of the strong lensing data we create 9 different 
reconstructions, each time removing one of the multiply 
imaged systems we use. The resulting A>:-maps do not 
change substantially, the main features (i.e. the elliptic- 
ity of the cluster and SW extension) seen in Fig.HJremain 
in all of the reconstructions. 

We have also run the reconstruction by changing the 
average redshift of weak lensing sources to zwl = 1.0 
and zwL = 2.0. The changes in mass estimates were at 
the 1% level; hence the precise value of zwl is unimpor- 
tant in this case. The normalisation of the mass profile 
is determined by the much stronger signal coming from 9 
multiple-image systems with "assumed" known redshifts. 
Weak lensing is helping here with the shape of the mass 
distribution in the areas where no multiply imaged sys- 
tems are present. 

Finally we study the depend ence of the re s ults o n ini- 
tial conditions. As noted in 'Bradac et al." (2006) the 
main features of the reconstruction are independent of 
the initial conditions we use. This is mostly attributed 
to the richness of the strong lensing data used in that 
paper and here. Still, we performed mass reconstruc- 
tions where we changed the velocity dispersion of the 
initial model to a = 1900 km s~^ and a = 1300 km s~^; 
the resulting mass estimates are within the errors quoted 
above. 

In summary, the errors on mass and surface mass den- 
sity quoted throughout this paper include the errors ob- 
tained from bootstrap resampling the weak lensing cata- 
logues, and removing individual strong lens systems (the 
latter dominate the error budget). The other errors dis- 
cussed in this section have a minor contribution to the 
total error budget. 

6.4. A Possible z = 4.08 System? 

During their spectroscopy run, Cohe n Kndbl ([200?) 
serendipitously discovered a candidate z = 4.08 object 
(see Fig. [6]). Using our final mass reconstruction we have 



searched for counter images of this object, unfortunately 
without success. If present we would have expected mul- 
tiple images to be easy to identify even in the absence of a 
reliable lens model, given its somewhat distinctive colors 
(see Figure [6]). Furthermore, the object is fairly bright in 
the F475W filter (mF475W,AB = 25.55 ± 0.03) and has a 
relatively blue F475W-F8i4W color (1.46 ± 0.04). This 
new piece of information appears to be inconsistent with 
the suggested redshift of z = 4.08, since the flux short- 
wards of Lyman al pha shou ld be absorbed by the Lyman 
alpha forest (Mad au et al.i ri996). 

At present we can therefore neither confirm nor reject 
the possibility of this object being lensed. However, we 
think it is unlikely that the bright knot is indeed at a red- 
shift of 4.08. Since we also see some underlying extended 
emission, it is possible that the bright knot is a fore- 
ground object and the extended object is being lensed; 
due to its low surface brightness we do not expect to be 
able to identify its counter images. A mo re likely possibil- 
i ty is that the emission line detected bv ICohen fc KneibI 
(|2002r ) is the [O II] doublet at 3727A, implying a redshift 
of z = 0.66 and therefore no multiple imaging. We do 
not expect the doublet to be resolved given the spectral 
resolution of ~ lOA(FWHM). We conclude that the red- 
shift of this object is uncertain at present and further 
investigations are needed to either confirm or reject the 
present redshift estimate. 

7. DISSECTING RX J1347.5-1145 INTO DARK MATTER 
AND BARYONS 

The relative distribution of dark matter and baryons 
in galaxy clusters is important for a number of reasons. 
Whereas baryons are only a minor total mass constituent 
in the dark matter dominated clusters, in the centres of 
the clusters baryons are an important mass component. 
In order to make a proper comparison with ACDM pre- 
dictions we need to disentangle both components. The 
interplay between baryons and dark matter is important 
and can potential ly change the inner slope of dark mat- 
ter halos (see e.g. iGnedin et al]l2004l ). Therefore in this 
section we separate the total mass measurements we ob- 
tain from gravitational lensing analysis into baryonic (gas 
and stars) and non-baryonic component and measure the 
dark matter distribution of RX J1347.5— 1145. 

7.1. 2-D distribution of dark matter and baryons 

To estimate the stellar mass distribution we first esti- 
mate the cluster K-band luminosity distribution by se- 
lecting the cluster members using the photometric red- 
shift (obtained using information from 8 colors from 
ground-based data in Bradac et al. 2005a) cuts at [0.3 — 
0.55]. We measure their K-band luminosity by as- 
suming the absolute solar r nagnitude Mk . ^t) = 3.28, 
galactic extin c tion f ollowing ISchlegel et al.l (|l998l ) and 
ICardel li et aT (1989") of Ak = 0.023, and K-correction 
from IPoggianti (1997) of KK(^d) = -0.25. We smooth 
the distribution using a Gaussian kernel of FWHM of 
30kpc. To convert the s t ellar l uminosity into stellar mass 
we follow iDrorv et al.l (|2004l ) and assume the stellar- 
mass-to- light ratio in K-b and to be M^/Lk = 0.74 ± 0.3 
(see also lBeh et alll2003[ ). 

The stellar K-band luminosity is shown in Fig. [H in 
white contours. It follows the total mass distribution 
considerately well, we detect the SW extension in both 



Fig. 4. — The F450W-F814W-F850LP color composite of the cluster RX J1347.5-1145. Overlaid in red contours is the surface mass 
density k, from the combined weak and strong lensing mass reconstruction. The contour levels are linearly spaced with Ak, = 0.2, starting 
at K = 0.5, for a fiducial source at a redshift of Zs oo. The X-ray brightness contours (also linearly spaced) are overlaid in yellow and 
the K-band light (tracing the stellar mass) is overlaid in white. North is up and East is left, the field is 2.3' x 2.3' , which corresponds to 
800 X 800 kpc^ at the redshift of the cluster. 



cases. Also shown is the X-ray surface brightness, in- 
dicating that the major baryonic component is spatiahy 
ahgned with the distribution of the total mass and hence 
dark matter within the uncertainties. 

7.2. Projected density profile of dark matter and baryons 

In Fig. [71 we present the total surface (projected) mass 
density profile S from the lensing reconstruction. Fur- 
ther we obtain the dark matter profile by subtracting 
the stellar and gas density profile from the total density 
profile, as described next. 

The stellar mass profile is obtained from the 2-D K- 
band light distribution. We assume the dominant source 
of error comes from the uncertainty of M^/Lk- We 
again smooth the distribution using a Gaussian kernel of 
FWHM of 30kpc; the exact value is of lesser importance 
here, since stars contribute a minor fraction of the total 



baryonic mass at radii > 50 kpc and we are not sensitive 
in smaller scale variations in the total mass profile. 

The projected gas profile is obtained from the observed 
X-ray surface brightness profile and the deprojected X- 
ray gas temperature profile. Finally, the dark matter pro- 
file is calculated from the difference between the baryonic 
and the total mass profile from strong and weak lensing 
mass reconstruction. In Fig. [71 we show the individual 
contributions of stars, gas, and dark matter. 

To estimate the inner slope of the dark matter halo we 
fit a generalized, projected NFW profile: the 3D density 
is given by 

/ \ PO,DM /-.x 

PBM[r) = ^ / ^o_g • 

(r/rs)^(l + r/rsr ^ 
We determine the asymptotic inner slope the scale ra- 
dius Ts and the normalisation po,DM from the data given 
in Fig. [71 We assume a flat prior > 0, since negative 
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Fig. 5. — The integrated, projected mass profile of 
RX J1347.5— 1145. The profile was determined by measuring the 
enclosed mass in cylinders, centered on the BCG using strong and 
weak lensing (hexagons), strong lensing using parametrised model 
(triangle) and X-ray data (crosses). We fit power law profiles to the 
strong and weak lensing data, the result is shown as a soli d line. 
1-sigma error region from the results obtained by Halkola ^et al.l 
(|2008j) using only strong lensing data is shown as dashed lines. 




Fig. 6. — 10^^ x 10^^ cut out of the z = 4.08 candidate discovered 
bv lCohen KneibI (|200l ) from the F450W-F814W-F850LP color 
composite. 

values of (3 give unphysical profiles. When fitting the 
profile we include the full covariance matrix (i.e. tak- 
ing into account correlations between the bins) for the 



estimate of the total density, gas density and the stellar 
mass density. The covariance matrix was determined di- 
rectly from strong and weak lensing bootstrap resampled 
reconstructions and by sampling the models given the un- 
certainties from gas and stellar mass measurements. We 
note however, that another potential source of system- 
atic errors will arise from the lack of the symmetry in 
the mass distribution and somewhat complicated geom- 
etry of this system. There is also the ambiguity when 
choosing the center of mass for the profile estimation, 
which should be computed around the center of mass. 
We chose to center the profile on the BCG, which is a 
good approximation in our case, since the maximum of 
S is consistent with the position of the BCG {2" ± 4^', 
i.e. < lOkpc). 

Furthermore, as already noted by ISand etal] (|2QQ7l ) 
there are strong degeneracies between the generalised 
NFW parameters. In particular we found a strong de- 
generacy between the scale radius rg and the inner slope 
p. The best fit profile gives P = 0.0 ± 0.1 and rg = 
160 kpc±10 kpc ( C200 = 15 and r2oo = 2400 kpc). How- 
ever, as shown in Fig. [71 the NFW profile fit with a fixed 
inner slope = 1 is not significantly worse. The latter 
gives a scale radius of = 350 kpcilOO kpc (c2oo = 6 
and r2oo = 2200 kpc). The statistical uncertainties alone 
on C200 are of the order of ~ 20%. 

To reliably determine the concentration parameter and 
hence /3 we therefore need the data to extend to much 
larger scales (i.e. beyond r^ir)- Strong lensing data alone 
(for this and also all the other clusters) are not sufficient 
for this task. We plan to extend this analysis by adding 
wide field weak lensing data that extend to a radius of 
~ 30' (10 Mpc); this will be a subject of a forthcoming 
paper. 

Finally, the thermodynamic X-ray pressure map 
(Fig. [3]) shows the localized region of high-pressure, 
shocked gas in the southeast quadrant. Therefore we are 
likely measuring a profile of two post-merging, unequal 
mass clusters. This is further supported by the presence 
of the eastern cD galaxy. Therefore the cluster might 
not be fully relaxed and therefore the measurement of 
the total dark matter distribution is unlikely to follow 
the NFW profile in great detail at these small radii. 



The longstanding puzzle of the discrepant mass esti- 
mates for the most X-ray luminous cluster seems to be 
definitely resolved. Both our strong and weak lensing as 
well as X-ray mass reconstructions agree well. In partic- 
ular we draw the following conclusions: 

1. Using the combined strong and weak lensing mass 
reconstruction we derive a high-resolution, abso- 
lutely calibrated mass map; we get projected, en- 
closed mass Ml(< 350kpc) = (5.9±0.5) x IO^^Mq. 
Within the same radius the projected mass de- 
rived from X-ray data alone is Mx(< 350 kpc) = 
(6.6t°;^) X IQI^Mq. 

2. The mass estimates are still in disagreement 
with previou s dynamical mass estimates by 
Cohen fc Kne ib (2002). It is however possible that 
Cohen fc Kneib (200 2) measured the velocity dis- 
persion of the infalling subcluster only () Allen et al.l 




CONCLUSIONS AND OUTLOOK 



11 




100 

R[kpcl 



Fig. 7. — (top) The projected mass density E profile of 
RX J 1347. 5— 1145 showing separately the stellar profile (filled tri- 
angles - green), the gas profile (stars - blue), dark matter profile 
(open triangles - black) and the total profile as measured from 
strong and weak lensing (hexagons - red). The dark matter profile 
has been fitted using the generalised NFW (cf. Eq. [T} - shown is 
the best fit model (f3 = 0, rs = 160kpc) and a best fit NFW model 
(f3 = 1^ rs = 350kpc). For gas and stars we use simple power law 
profiles. Also shown for reference are the E-profiles for 3D density 
profiles p (X and p oc r~^. Dark matter points have been offset 
in the R-direction for clarity, (bottom) The ratio of dark matter 
to total matter. 

|2QQ2[ ), thereby biasing the mass estimate low. If 
taken as a measurement of the cluster mass, such 
a low velocity dispersion cluster is in complete dis- 
agreement with strong lensing features we observe 
(regardless of their redshift). We are currently 
analysing the remainder of spectroscopy data to 
resolve this issue (Lombardi et al. in preparation). 

3. Using the exquisite resolution of ACS data we are 
able to study the spatial distribution of dark mat- 
ter with respect to the baryons at unprecedented 
accuracy. We clearly detect a mass concentration 
centered on the BCG and a SW extension, which 
follows the light distribution of the cluster mem- 
bers. Further we fit a generalised NFW model to 



the dark matter density profile, finding strong de- 
generacies between the inner slope /3 and the scale 
radius r^. Our data does not extend to sufficiently 
large radii, for more secure determination of rg, 
and concentration parameter; we plan to use weak 
lensing data extending beyond the virial radius in 
the future to resolve this issue. 



The cluster RX J1347.5— 1145 is an extremely valu- 
able object for understanding the details of cluster for- 
mation and evolution. In addition, due to their magni- 
fying power, these most massive clusters are ideal tools 
to study the high-redshift universe. By searching around 
the critical line for high redshift sources as predicted by 
our lens model we were already able to find high redshift 
multiple imaged systems. We plan to extend this search 
in the future. 

The large wavelength coverage of the data for this clus- 
ter has proved to be extremely valuable for a detailed 
study of the cluster mass distribution. With such data 
sets most massive clusters of galaxies remain one of the 
key objects to constrain cosmological parameters and to 
study the galaxy formation and evolution from the early 
times until the present. 
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